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1.  SYSTEM  MODEL 


The  front  end  of  the  radar  receiver  is  shown  below: 


Let  the  voltage  transfer  function  relating  the  IF  amplifier  output 
signal  to  the  RF  amplifier  input  signal  be  the  zero-memory 
nonlinearity  specified  by 

21  f*  -Z2/2CT2 

v  =  g  (v.)  -  -==.  e  g  dZ . 

01  <2^  J 

9  0  . 


Due  to  the  mixer,  the  carrier  frequency  of  vq  is  translated  from  that 
of  v^  by  an  amount  equal  to  the  local  oscillator  frequency.  Other 
than  for  this  frequency  shift,  the  spectrum  and  correlation  function  of 
vq  is  unaltered.  For  simplicity,  the  nonlinearity  is  modeled  as 
though  vq  and  v^  have  identical  carrier  frequencies.  The  voltage 
transfer  function  accounts  for  the  cumulative  nonlinearities  in  the  RF 
amplifier,  mixer,  and  IF  amplifier. 
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Observe  that 
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Hence,  the  voltage  transfer  function  is  modeled  as  a  soft  nonlinearity 
with  saturation  levels  at  +t.  Note  that  the  integrand  of  g(-)  is 
proportional  to  a  zero-mean  Gaussian  probability  density  function 
(pdf)  with  variance  o^.  is  referred  to  as  the  rms  value  of  the 

voltage  transfer  function.  Let  a i  denote  the  rms  value  of  v^  The 
ratio  of  these  two  rms  values  yields  the  parameter 


When  v.  -  a  ,  the  output  is 

i  g 


>" g  “V  ‘  vs?  F 
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e  y  dy  -  0 . 68t .  (5) 


The  slope  of  the  voltage  transfer  function  is 


o  _  2t 

dvi  "  V2Jtor 
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l2  t 

At  the  origin,  therefore,  the  slope  is  yjf  Consider  the  linear 

transfer  function. 
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V  -  f  (V,  ) 
O  1 


(7) 


When  v.  ■  a  ,  the  output  of  the  linear  transfer  function  is 
i  g 


fi 

-vji  t  »  0.80  t. 


(8) 


Hence,  when  v.  *  <T  ,  the  nonlinear  transfer  function  is  down  from  an 
i  g 

ideal  linear  transfer  function  by 
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.  68£ 

10  .soe 


1.4  dB. 


(9) 


In  this  sense,  is  a  measure  of  the  linearity  of  the  voltage 
transfer  function.  On  the  other  hand,  CT  is  a  measure  of  the  input 
fluctuations.  Thus,  a  is  an  indicator  of  the  "width"  of  the  linear 
portion  of  the  voltage  transfer  function  relative  to  the  rms  value  of 
the  input  signal. 


2.  DISTORTION  OF  THE  PDF  CAUSED  BY  THE  NONLINEARITY 


Let  the  pdf's  of  and  Vq  be  denoted  by  pv  (v^)  and  pv  (v  ) , 

i  o 


respectively.  It  follows  that 
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From  Page  2, 


Vi"g"1(Vo)- 


(10) 


dv.  V2  no2  +v*/2  a2 
_±  _ _ 2_  e  1  g 

dv  2 1 

o 


(11) 


Let  be 


a  zero-mean  Gaussian  random  variable  with 


variance  <7 
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Then 


PVi<Vi> 
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It  follows  that 
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Simplifying,  we  have 
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For  a  -  1, 


Pv  <V 
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0  ;  otherwise. 


oCz) 


Thus,  the  IF  amplifier  output  is  uniformly  distributed  when  a,  ■  o  . 

For  a  >1,  the  exponent  in  Pv  (v  )  is  negative.  Note  that 

o 


g  1(f)  -  00  and  g  ^(-t)  -  . 
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Thua,  the  density  function  goes  to  zero  at  vq  -  +£.  The  output 

density  function  is  approximately  Gaussian  provided  a  »1.  However, 

irrespective  of  how  large  is  a,  p„  (v  )  can  never  be  Gaussian  since  v 

Vo  0  ° 

is  constrained  to  the  interval  (-£,£)  whereas  a  Gaussian  random 

variable  can  assume  values  in  the  infinite  interval  (-<»,  °°)  .  For 

a  <  1,  the  exponent  in  p^  (v  )  is  positive.  The  density  function  then 

o 

goes  to  <»  for  vq  =  +  £ .  In  the  limit,  as  a  -»  o,  the  density  function 
approaches  two  impulses  as  shown  below. 


0  and 


(19) 


At  the  origin,  g  (0) 


PV  <V  ■  21 


3.  DISTORTION  OF  THE  CORRELATION  FUNCTION  CAUSED  BY  THE  NONLINEARITY 
Define  the  Fourier  transform  of  the  voltage  transfer  function  to 


OO 

-juvi 

ju)  -  g(v^)  e  dv^ 


The  output  correlation  function  is  given  by 

R  (X )  -  E[v  (t)  v  (t-X)]  . 

o  o  o 

The  output  voltages  can  be  expressed  in  terms  of  the  Fourier 
transform  of  the  voltage  transfer  function  according  to 
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It  follows  that 
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The  quantity 


M( ju,  jw) 
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f  juvi(t) 

e  e 


jwv^  (t-X ) 


(26) 


is  recognized  as  the  joint  characteristic  function  of  the  random 
variables  v^(t)  and  v^(t-T).  Hence, 


R  (X)  -  — 2 
O  4JI 


G(ju)  G(jw)  M ( ju,  jw)  du  dw.  (27) 


— oo  —  oo 


Let  v^(t)  and  vMt-X)  be  jointly  Gaussian  zero-mean  random 

variables  with  variance  a2  and  correlation  function 

9 

(X )  -  E  (v^t)  vi(t-T)  )  . 

Their  joint  characteristic  function  is  then  given  by 


(28) 


M ( ju,  jw)  -  e 


where 


-1/2  [<T;  u2  +  2R.  (X)  u  w  +  O2  w2  ] 
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R.  (X) 

Pi(X)  -  ~P' 


(29) 


(30) 


(31) 


is  the  correlation  coefficient  relating  v^(t)  and  v^(t-T).  As  a 
result,  the  output  correlation  function  can  be  expressed  as 


R  (x ) 
o 


OO  OO 

«!  I  IG,3u 


-«T;  / 2)  [u2  +  w2  ]-R,  (X  )  uw 
)  G ( jw)  e  1  1  dudw.  (32) 


To  evaluate  G(ju),  we  make  use  of  the  differentiation  property  of  the 
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Fourier  transform.  To  begin  with,  we  have  the  Fourier  transform  pair 


g(vi) 
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V2j to2 


-z2 /2a2 
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Using  the  differentiation  property, 
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e  i  g 


ju  G ( ju) . 


However,  it  is  well  known  that 


-v2  /2a2 
e  i  g 
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It  follows  that 
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Thus, 
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Substituting  into  the  expression  for  Rq  (X )  ,  we  have 


*  «>  - 
o  4Jt 


-«J2/2)[u2  +  w2]  -  (CT2. 12)  [u2  +  w2]-R,(t)uw 
(l/uw)e  ^  e  1  dudw. 


Combining  terms  in  the  exponent  of  the  integrand. 
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To  obtain  an  integrand  without  the  factor  of  —  ,  we  differentiate 

uw 


R  (X )  with  respect  to  R. (X) .  Specifically, 
o  r 
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Now  make  the  change  of  variables 


la2  +CT2 


u  and  y 


a2  +CT2 
_2 _ i 


(41) 


:hi3  results  in 
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To  further 'simplify  the  intregand,  introduce  the  variables  V  and  7 
such  that 
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v  =  x  +  ay 
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where 


Ri(T) 

3  " 
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The  Jacobian  of  this  transformation  is  given  by 


(43) 


(44) 


(45) 


3h 

d  h 
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3v 
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- 
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Hence, 


dxdy  *  7:  — ,  dv  dY 

Vl-a2 


(47) 


With  this  transformation,  the  exponent  of  the  integrand  becomes 
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The  double  integrand  is  readily  evaluated  to  yield 


JJ-- 


-(v2+y2) 
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(52) 


Consequently, 
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Recalling  that 


d  (sin  1  u) 
dx 
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du 
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it  follows  that 
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cr2  +02 
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where  c  is  a  constant  of  integration.  To  evaluate  c,  consider  the 
special  case  for  which 
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and  X  m  0 .  Then 
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1.  From  page  4  of  the  notes,  the  output  pdf  is 
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Hence, 
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It  follows  that  c  ■  o.  The  output  correlation  function,  therefore,  is 


given  by 


(21) 2  .  -1  Ri(X) 
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Recalling  that 
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and  a 


R  (t)  can  also  be  expressed  as 
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. 

4 .  DISTORTION  OF  THE  POWER  SPECTRAL  DENSITY  CAUSED  BY  THE  NONLINEARITY 


The  power  spectral  density  at  the  output  of  the  nonlinearity  is 


II 


given  by 


S  «D) 
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R  (X)  e  dT 
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i2tj  i  r  .  -1  pi 1,1  -j“ 

~2T  sin  TS3-  e  <“• 


Assume  the  normalized  input  correlation  '  .  ction  is 


-BIT  | 


Pi(T)  -  e 


This  corresponds  to  white  Gaussian  noise  passed  through  a  lowpass  RC 
filter  of  bandwidth  B.  Hence,  the  power  spectral  density  at  the 


nonlinear  input  is 


„  ,  .  2B  2  1 

Si(W)  "  B2  +CJ2  ”  B  1+  (CO/B)  ‘ 


To  evaluate  the  output  power  spectral  density,  we  expand 

-1  i 

sin  [pi (t ) / (1+a  ) ]  into  a  Taylor  series  and  integrate  term  by 
term.  For  simplicity,  let  a  -  1. 


Then 
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By  inspection,  the  output  power  spectral  density  is  given  by  the 
series  expansion 
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iMii  i 

2K  B 


1+  {(D/B) 2 


+  .014 


1+  (Q)/3B) 


T  +  .000936 


1+  (CO/ SB) 2 


(73) 


In  general,  the  output  spectrum  can  be  thought  of  as  resulting  from 
the  superposition  of  an  infinite  number  of  spectra  having  bandwidths 
B,  3B,  5B,  etc.,  as  shown  below. 


Because  the  infinite  series  converages  rapidly,  the  total  power 
outside  the  interval  (-B,  B)  is  relatively  small.  Nevertheless,  it 
may  be  significant  in  the  subclutter  visibility  application. 
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Abstract 

This  paper  deals  with  the  compensation  of 
mutual  coupling  effects  in  a  sensor  array  vhen 
using  the  subspace  approach  known  as  the  matrix 
pencil.  This  method  is  easy  to  implement  and  is 
applicable  even  in  the  presence  of  coherent 
sources.  Computer  simulation  shows  that  use  of  the 
mutual  Impedance  matrix  effectively  compensates 
for  the  presence  of  mutuals. 

Introduction 

Direction  finding,  which  involves  the 
estimation  of  the  angles  of  arrival  of  sources,  is 
very  important  in  many  sensor  systems  such  as 
radar,  sonar,  seismology,  etc.  Over  the  years 
several  methods  have  been  proposed  to  solve  this 
kind  of  problem  [1,2].  Recently,  subspace  ap¬ 
proaches  have  been  introduced  (2-7|.  They  are 
based  on  an  eigenvalue-eigenvector  decomposition 
of  the  correlation  matrix.  They  are  shown  to  yield 
asymptotically  unbiased  estimates  of  such  signal 
parameters  as  angles  of  arrival  and  number  of  sig¬ 
nals.  However,  these  methods  have  not  taken  into 
account  effects  of  mutual  coupling  between  array 
elements  which  can  significantly  alter  the 
eigensystems  underlying  the  solution  procedures. 
In  this  paper,  we  are  dealing  with  the  compensa¬ 
tion  for  mutual  coupling  effects  when  using  the 
subspace  approach  known  as  the  matrix  pencil  ap¬ 
proach.  Computer  simulations  show  that  significant 
improvement  in  the  estimates  for  the  angles  of  ar¬ 
rival  can  be  achieved. 


nj(t)  is  the  additive  noise  assumed  to  be 
zero-mean, 

*  wDsinO^l/c 

o i  is  th  center  frequency  of  each  of  the  spa¬ 
tial  sources. 

c  is  the  speed  of  propagation  of  the  plane 

waves , 

a(9)  is  the  beam  pattern  of  each  sensor. 


Given  the  number  of  sources  d  and  the  m  data 
points  ,  y i ( t , 0)  ,  we  form  (d+1)  vectors  Yn  of 
length  (m-d),  where 

Xn^  *  (yn’^n*! '  •  •  •  •  »yn*m-d-l)  '<  •  •. (<•♦!). 


The  two  matrices  and  Nj  are  then  formed 
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'  T  T 

1  1 

t 

1 

M1  * 

Xl  12  • 

•  Xd 

i  Nj  - 

X2  X3  • 
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1  1 

L  4  A 
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A 

1  1 
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1 

A 

Let 


.  .  .  1 

.  .  .  e^d 


ejfm-d-Ddj  ej(m-d-l)*2 


ej(m-d-l)*d 


Moving  Window. 

Assume  we  have  a  linear  array  composed  of  m 
Identical  sensors  with  uniform  spacing  D.  Let 
there  be  d<  m  narrowband  sources  located  at 
azimuthal  angles  6^  i  k«l,2,...,d  ,  emitting  sig¬ 
nals  whose  complex  envelopes  are  denoted  by  sk(t); 

k«l,2  ,3 . d.  Assume  the  sources  are  in  the  far 

field  such  that  planar  wavefronts  arrive  at  the 
array.  With  reference  to  the  first  sensor,  the 
received  signal  at  the  i**1  sensor  is  modeled  as 

d 

y.(t,e)-E  sk(t)a(9k>«j<i'1)*k  *nj(t)  (1) 

~  k.l 

for  i-I . m, 


where 
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ST  -  lsltS2, - ) , 

and 

t*n"^"  *  I  nn,nn*  1  >  •  •  •  •  »nn*m-d- 1  I’ 
Then  Yn  can  be  rewritten  as 
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xn  -  A4 

BS 
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E|H1llM1)=UliVU 

♦(m-d)a2  Id 
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matrices 

and  Nj 
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where  Id  is  the  (dxd)  identity 
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Let  F,  S,  U,  N' ,  and  n"  be  the  matrices 
t  t  t 

I  I  I 

S  BS  .  .  . 

I I  I 

1  i  .  4 


u  - 


1  eJ*l  . 
1  eJ*2  . 


sd 

.  el**'1)*! 


l  i  eJ*d  .  .  .  eJM-D* 


T  T 


i!i  52 


5d 


I  I  I 

14.  1 


(4) 

(5) 
'id 


Defining 

Fpq  -V  eJ(i-lHV*q>. 
i-1 

Spq  -  E|sq*P  J' 


the  matrix  V  can  be  vritten  as 


sllaUFll  s21a21F21 
s12a12F12  s22a22F22 


sldaldFld  s2da2dF2d 


sdladlFdl 

sd2ad2Fd2 


sddaddFdd 


Thus,  if 

M  .  E|djHMj]-(m-d)o^  Id  -UgVU  (6) 

N  -  ElN^HxMm-d)*2  lld-Ug*"vu.  (7) 

The  matrix  pencil  then  becomes 

H-XN  -  UgVU-XUgB^VU  -  UH(I-X*H)VU  (8) 

which  satisfies  the  requirements  of  the  pencil 
theorem.  Hence,  the  values  of  X  for  which  the  rank 
of  M-XN  decreases  by  1  are  given  by 


Xu  »  el^k  ;  k«l,2, . . .  ,d. 


(9) 


t  t 

I  I 

52  53 


1  1 


5d*l 


Then 


.  A8SU*N  and  N1  •  ABSBLUN 


(3) 


Assuming  the  signals  and  noise  to  be  statistically 
Independent  and  that  the  noise  components  are  un¬ 
correlated  from  sensor  to  sensor,  we  get 


The  angles  of  arrival  are  given  by 

©jt  »  sin'^-jclnfXj^/uO) ;  k«l,2,...,d.  (10) 

Mutual  Coupling 

Consider  a  linear  array  of  m  dipoles 
uniformly  spaced  at  a  distance  D.  Each  dipole  is 
of  length  t  and  has  a  radius  r  satisfying  the  con¬ 
dition  r<<t.  A  load  is  attached  to  the  center  gap 
of  each  dipole.  Assume  there  are  d  narrowband  sig¬ 
nals  impinging  on  the  array  as  planar  wavefronts. 
The  voltages  induced  by  the  assumed  signals  on  the 
loads  are  the  outputs  of  the  dipoles.  Induced  cur¬ 
rents  will  appear  on  the  dipoles.  These  currents 
reradiate  and  generate  scattered  fields.  The  scat- 
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In  matrix  notation 


tered  fields  then  induce  currents  on  the  neighbor¬ 
ing  dipoles.  The  process  of  Induction  and  reradia¬ 
tion  causes  the  mutual  coupling  among  the  dipoles. 

Using  one  sinusoidal  expansion  and  weighting 
function  per  dipole,  the  method  of  moments  (9,10) 
was  used  to  obtain  the  matrix  of  mutuals.  Denote 
the  current  distribution  by  J(z)  (assuming  longi¬ 
tudinal  distribution  and  neglecting  all  other  dis¬ 
tributions)  and  the  j-th  expansion  function  by 
ft(t). 

Then 


V<s>=-Z  X  (20) 

where 

y(s)T,|V(s)(1)iV(s)(2) . V<s)(m)  j 

and 

IT=[ I ( 1 ) , I ( 2) , .  .  • , I (m) ] . 

The  matrix  Z  can  be  decomposed  into  two  parts  as 


m 

J(z).I  KjjMz)  (11) 

1-1 


where  I(J)  are  unknown  amplitudes  to  be 
determined.  At  a  point  (y,z)  in  the  '/-Z  plane,  the 
scattered  field  is  given  by 

m 

E(s>(y,z).r  I(j)E<5|y,z)  (12) 

1-1 

where  Ej(s)(y,z)  is  the  scattered  field  from  the 
j-th  dipole.  The  total  field  will  then  be 


where 

and 


Z=Zq+Zl 

Zq  is  the  generalized  impedance  matrix 


Z(_  is  the  load  matrix. 

Assuming  that  all  loads  are  loaded  with  the  same 
load  z^,  the  matrix  ZL  is  given  by 


ZL* 


0 


E(y,z)-B(inc)(y,z)  ♦  E(s>(y,z)  (13) 

where  E^nc^  is  the  incident  field.  Let  Ez  be  the 

z-component  of  the  total  field.  A  generalized 
voltage  V(i)  induced  on  the  subsection  spanned  by 
the  function  fj(z)  can  be  defined  with  respect  to 
a  weighting  function  w(z)  as 


L  zi  1 

The  Ij-th  element  of  Z,  therefore  is 

ziJ-ZiJ«l»ij- 

The  voltages  induced  on  a  load  z^  are  given  by 


V(i)-F(Ez(y,z),wi(z))  (14) 

where  F  is  bilinear  with  respect  to  Ez(y,z)  and 

"{(y.z)- 

Similarly,  we  define 

v(inc>(i)=F(E£inc)(y,z),wi(z))  (15) 

and 

V<s)(i)»F(Eis>(y.z),w1(z)).  (16) 

Thus, 

V(i)=.v(lnc>(i)*v(s)(l),  which  becomes,  for 
metallic  scatterers, 

V(l).V<lnc>(i)*V<s>(l).0 

or 

V(lnc)(i),_v(s)(i)  (17) 

However, 

m 

V<s>(i)  .F(  E  I(j)E<sJy,z),wi(z)) 

J-l 

m 

-E  Z(j)F(E<fiy,z),wi(z)) , 

j-l 

Let  z*j  be 

i1J.-F(E<jiy,z),vi(z)).  (18) 

Thus 

m 

V<*>(1).E  -z 1 j  I(j)  !  1.1,2 . m.  (19) 

J*1 


V(f>  =  ZL  I  and  I  -  ZL_1  V(t). 


However , 


v<  .  ZI  .  Z0ZL ♦  V(t), 


which  implies  that 

v< 1 >  - ( i *z0zL_1 l _1  v(inc).  (21) 

Let  H  be  the  matrix 


H=  IIvZqZl-1]- 
H  can  be  writ  ten  as 


(22) 


U(zu/Zl)  (z12/ Zi)  ... 

<z21/zl>  U(z22/z1)  .  .  .  (z2m/zi) 


(zml/zl>  <zm2/zl> 


1  *■!  zmm/zl  1 


« 


Thus,  when  incident  signals  are  impinging  on  the 
array  and  in  the  presence  of  additive  noise,  the 
output  of  the  linear  array  will  be 

v(D  =H"  1  V(inc>  *  N. 


For  simplicity,  let 
x,y( inc) 

and 

Y=y( r > . 

We  now  have  a  relationship  between  the  incident 
signals  and  the  received  signals  at  the  outputs  of 
the  array,  which  is 


Y  .  H'1  X  ♦  N. 


(23) 
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If  we  try  to  use  the  vector  Y  in  the  formulation 
of  the  pencil  theorem,  it  is  not  possible  to  ob¬ 
tain  the  decomposition  needed  in  that  formulation. 
An  estimate  X  of  X  is  needed  in  order  to  obtain 
the  decompositions  needed  in  formulation  of  the 
pencil  theorem. 

Assuming  that  the  signals  and  noise  are 
statistically  independent  and  that  the  noise  com¬ 
ponents  are  uncorrelated  zero-mean  random  vari¬ 
ables  vith  variance  e  ,  the  minimum  error  mean- 
squpred  linear  estimator  results  when  the  error 
(X-X)„  is  orthogonal  to  the  observed  data  Y. 

Let  X>  R  Y,  where  R  is  to  be  determined-!  After 
some "compu tat  ion  and  defining  C  to  be  the  correla¬ 
tion  matrix  of  the  observed  data  Y,  we  obtain 

R  -  H  (C-o2  I„)  C'1.  (24) 

Computer  Simulation 

The  scenario  used  for  this  simulation  consisted  of 
two  coherent  sources  (d»2)  which  are  incident  on  a 
linear  array  consisting  of  eight  half  wavelength 
dipoles  (m>0).  The  sources  are  assumed  to  be  lo¬ 
cated  at  8j«18*  and  02*22°  .  The  noise  was  simu¬ 
lated  to  be  vhite  Gaussian  vith  zero-mean  and  unit 
variance.  The  sensors  vere  positioned  at  half 
vavelength  apart  such  that  uD/c  *  n  .  The  results 
of  the  simulation  are  shown  belov. 


(Vithout  compensation  for  the  mutuals) 


SNR 

mean 

9i 

mean 

02 

variance 

9i 

variance 

«2 

30  dB 

20.6537 

41.4553 

4.177E-2 

1.98405 

25  dB 

20.6478 

41.5294 

4.556E-2 

2.98911 

20  dB 

20.6374 

41.6539 

5.787E-2 

7.64255 

15  dB 

20.6205 

41.1151 

1.192E-1 

44.5549 

(  100  snapshots/run,  50  runs) 

(With  compensation  for  the  mutuals) 

SNR 

mean 

9i 

mean 

e2 

variance 

6i 

variance 

e2 

30  dB 

17.9308 

21.9435 

5.164E-3 

8. 784E-3 

25  dB 

17.8171 

21.7513 

1.947  E-2 |  3.909E-2 

20  dB 

17.5739 

20.9331 

1.486  E-l 

2.019E-1 

15  dB 

14.1798 

20.1921 

15.63312 

5.963E-2 

(  100  snapshots/run,  50  runs) 


Note  that  extremely  poor  estimates  are  obtained 
without  compensation  for  the  mutuals.  Compensation 
results  in  significant  improvement. 
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Abstract 


Problem  Formulation 


The  moving  vindov  (1|  is  one  of  several 
operators  that  can  be  used  in  conjunction  vith 
Che  matrix  pencil  approach  for  estimating  the 
locations  of  multiple  sources.  A  major  advantage 
of  the  moving  vindov  is  that  the  locations  of  d 
(d<a)  completely  correlated  narrowband  sources 
can  be  estimateo  using  a  linear  equispaced  array 
of  m  sensors.  In  this  paper,  the  concept  of  the 
moving  vindov  is  extended  to  the  wideband  case 
by  using  the  Fourier  series  approach  to  convert 
each  wideband  signal  into  a  sum  of  narrovband 
signals.  A  transformation  matrix  is  then  intro¬ 
duced  Co  generate  a  single  equivalent  narrovband 
signal.  Any  narrovband  high  resolution  technique 
can  be  used  at  this  point  to  generate  the 
estimates.  Computer  simulations  show  that  the 
muring  vindov  performs  better  than  ESPRIT. 

Introduction 

Many  high  resolution  algorithms  have  been 
developed  for  the  case  of  narrovband  signals. 
Hovever,  narrovband  modeling  is  not  appropriate 
for  videband  signals.  In  particular,  the  time 
delay  cannot  be  approximated  by  a  phase  shift. 
Recently,  several  veil  known  narrovband  high 
resolution  algorithms  have  been  extended  to  the 
videband  case  (2-6|. 

An  effective  approach  consists  of 
decomposing  the  videband  signals  into  sums  of 
narrovband  signals  l 3—5 ] .  Wax  et.  al.  [3]  use 
the  HUSIC  algorithm  to  obtain  a  spectral 
estimate  by  averaging  either  geometrically  or 
arithmetically  the  contributions  from  each  of 
the  subbands.  This  incoherent  processing  is 
referred  to  as  post  averaging.  Wang  and  Kaveh 
(4 |  have  proposed  a  pre-processing  scheme  using 
linear  transformations  to  combine  the  various 
sub-bands  into  a  single  band.  As  in  the  previous 
paper,  the  HUSIC  algorithm  is  employed.  However, 
nov  only  one  eigendecomposltion  is  needed  to 
estimate  the  angles  of  arrival.  This  approach 
knovn  as  coherent  signal  subspace  processing 
(CSS)  is  shown  to  outperform  the  previous  one. 
Hovever,  a  disadvantage  is  that  a  preliminary 
crude  estimate  of  the  emitter  locations  is  re¬ 
quired.  If  the  angles  are  clustered  within  a 
beamvldth,  the  method  performs  well.  Otherwise, 
spatial  prefiltering  is  needed.  The  contribution 
of  this  paper  is  to  shov  that  CSS  can  be  com¬ 
bined  vith  the  moving  window  operator  [1|  to 
solve  for  the  angles  of  arrival.  Computer 
simulations  reveal  that  the  newly  proposed  tech¬ 
nique  performs  better  than  ESPRIT  (7). 


Consider  the  problem  of  estimating  the 
angles  of  arrival  of  wideband  signals.  The  no¬ 
tion  of  Fourier  coefficients  is  used  here  in 
conjunction  with  the  matrix  pencil  approach.  As¬ 
sume  that  all  incoming  signals  have  approximate¬ 
ly  the  same  bandwidth  B.  Let  T  be  an  observation 
interval  and  let  ano  f^  be  the  lowest  and 
highest  frequencies  contained  in  B.  In  practice, 
this  band  is  determined  by  Fourier  decomposition 
of  the  received  signals.  Assume  ve  have  a  linear 
array  composed  of  m  identical  wideband  sensors 
uniformly  spaced  at  a  distance  A.  Let  there  be  d 
(d<m)  videband  sources  located  in  the  far  field 
so  that  plane  waves  arrive  at  the  array.  The 
received  signal  at  the  i-th  sensor  can  be 
modeled  as 

d 

x1(t)-Ea(ek)sk(t-tlk)+ni(t)}  i-1 . .  (1) 

k-1 


where  Tjk  is  the  time  delay  that  source  k  takes 
to  travel  from  the  reference  point  to  the  i-th 
sensor.  Taking  the  reference  as  the  first 
sensor,  Tjk  can  be  written  as 

Tlk.(l-l)(4/c)sin(^)  (2) 

where  c  is  the  speed  of  propagation  of  the 
waves.  Define  the  Fourier  coefficients  of  the 
signal  received  at  the  i-th  sensor  by 


Xt(ur)=  (T)v» 


T/2 

xt(  t)exp(-jurt)  dt 
-T/2 


(3) 


where 

wr-(2n/T)(r1*r)  ,  r=l,2,.  .  .,R, 

rj  is  a  constant  and  R  is  the  number  of  sub¬ 
bands. 

With  reference  to  equation  (1),  the  r-th 
Fourier  coefficient  of  xj(t)  can  be  expressed  as 

d 

Xi(»r)-!:a(ek)Sk(u)r)e-j(l-1),t,k(“r)+N1(tor)  (4) 

k-1 

for  1  a  1 , 2 , .  .  .,m, 

where  the  r-th  Fourier  coefficients  of  sk(t)  and 
nj(t)  are  denoted  by  Sk(wr)  and  Nj(»r),  respec¬ 
tively,  and 

*k(«r>— KHA/cJsinOic).  (5) 
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(9) 


( d 1 )  vectors  Xn(ur)  o£  length  <m-d)  are  formed 

where 

•  •  •  *n+m-d-l(“r) • 
n-1,  2  . . .  ,(d+l) 

It  Is  easy  to  show  that  X,,^)  can  be  put  in  the 
for* 

Xn(Mr)-A(»r)»<n_1)(«r)BS(«r)  ♦£)„(&>£.),  (6) 

whera  A(ut),  *(«r),  B,  §(«,.)  and  N^Wf)  are 

1  ...  1 
eJ*l<“r)  ...  e-l*d<“r> 

Am 

^(■-d-l)*!^)  ...  eJ(m-d-l)*d(wr) 


R 

Xn(ooo)  =  (1/R)  Z  Tn(ur)Xn(«r). 
r  =  l 

Let  S'  and  N^,  be 
R 

S'  -  (1/R)  l  Tn(ur)S(wr)  ,  (10) 

r»l 


R 

-  (1/R)  Z  Tn(wr)Nn(or)  .  (11) 

r*  1 

Therefore,  Xn(<^j)  can  be  expressed  as 

Xn(rJ0)=A(U0)*(n-1>(u0)BS'-N^.  (12) 

In  the  remainder  of  this  paper  the  dependence  on 
uq  is  assumed.  The  two  matrices  Mj  and  N|  are 
then  formed  where 


eJ*l<“r> 

eJ+2<“r)  0 

*  . 

°  ’ej*d<“r> 

ST  -  tS1(«r)  S2(o)r)  .  .  .  Sd(ur)  I, 
and 

SJn"^  ”  lNn(“r)  •  '  '  ^n+m-d-l^ur^  I  * 

Transformation  matrices  Tn(wr)  (3|  are  then  used 
in  such  a  way  that 

Tn(Ur)A(tor)*<n-1>(ur)-A(U0)*(n-l)(o)0)  (7) 

where  uq  is  a  conveniently  chosen  frequency, 
usually  selected  as  being  the  midband  frequency. 
Then 

Tn(Ur)Xn(Ur).Tn(Wr)A(«r)*(n-1)(wr)8S(u.l.) 

*  Tn(^)Nn(»r) 

=A((a0)*<n-1)(«0)BS(»r)*Tn(ur)Nn(«r).  (8) 

Assuming  that  we  have  R  sub-bands,  we  form  the 
arithmetic  average  as  given  by 
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1  e-Xl  .  .  .  « J  (d - 1 ) ♦! 
1  e-)*2  .  .  .  e-)(d-l)*2 


i  *J*d  .  .  .  ^(<*-i)*d. 


N 


t  t  T 

i  I  I 

Si'  32'  •  •  •  3d' 

!  I  I 

4  4  A 


r  t 


T 


N 


32'  33'  •  • 

I  I 

4  4 


Then 

Hj  -  ABSU+N' 
and 

-  ABSKI+n". 


3d*l' 


(15) 

(16) 


Assuming  that  the  signals  and  noise  are 
statistically  independent  and  that  the  noise 
components  are  uncorrelated  from  sensor  to 
sensor,  we  get 


ElM^MJ-U^U  ♦  E[N'hN'|  (17) 

E[N1HMi|.UH*HVU  +  E|N"HN' |.  (18) 

where  V  is  the  matrix  V»E[SHBHAHABS| .  Defining 


m-d 


eJd-lXVV. 


rite  matrix  pencil  then  becomes 

M-XfJ  =  UHVU-XUH*HVU  =  U,,(I-X*H)VU  ( 21 ) 

which  satisfies  the  requirements  of  the  pencil 
theorem.  Hence,  the  values  of  X  for  which  the 
rank  of  M-XN  decreases  by  1  are  given  by 

xk  =  eJ*k  ;  k=l ,  2 . d.  (22) 

The  angles  of  arrival  are  given  by 

®lt  =  sin*^{jcln(Xk)/ojQO)  ;  k-1,7 . il  (21) 


Simulation 

Several  possibilities  exist  for  choosing 
the  transformation  matrices  Tn  |7|.  It  can  be 
shown  that  a  diagonal  transformation  leads  to 
the  simplest  analysis.  Assuming  the  sources  to 
be  clustered  within  the  proximity  of  one  loca¬ 
tion  6.  the  transformation  matrices  Tn(wr)  then 
become 

Tn(«r)=e-j(n-l>(«o-"r><A/c)sin(3)  Tj(«r) 
where 

T ! ( wr )  -  diag{  1  e_j (<^o-“r )( ^c)sin(0) 

.  .  .  e- j  (m-d-1  )(u)g-ur )( A/c)sin(B)  j  _ 

With  this  transformation  it  follows  that 

Tn(»r)A(wr)*<n~l>(wr)  =  A(«0)*(n-1>(fctt). 

Assuming  that  the  noise  components  are  uncorre¬ 
lated  from  sensor  to  sensor  and  from  sub-band  to 
sub-band  with  tern  mean  and  variance  o’  ,  i t  can 
be  shown  that 

E|N'hN' |  =  R  o2  Id 

E1N"hN'  ]  =  R  o2  I]d 


Spq  -  E[S*Sp  1, 

apq  *  aqV 

the  matrix  V  can  be  written  as 


sllallFll  ....  SdladlFdl 
s12a12F12  ....  Sd2ad2Fd2 


sldaldFld  ••••  SddaddFdd 


Define  the  matrices  H  and  N  as 


H  -  E|M1HM11-EIN'HN' ]  -  UH  V  U 

and 

N  -  E(NjHM1|-E(N"hN'  ]  -  UH  V  U. 


(19) 

(20) 


where  Id  is  the  dxd  identity  matrix  and  Ijd  is 
the  matrix 


0  10  0 
0  0  10 
0  0  0  1 


lid  - 


0  0  0  0 
0  0  0  0 


0 

0 

0 


1 

0 


In  the  simulations,  we  have  considered  a  linear 
array  consisting  of  8  sensors  uniformly  spaced 
at  a  distance  6*c/(2  fQ>  .  Following  the  example 
in  (7),  the  two  sources  were  assumed  to  be  lo¬ 
cated  at  16’  and  20’  and  to  have  ideal  rec¬ 
tangular  spectra  of  bandwidth  B.40  Hz  centered 
at  f0-100  Hz.  The  broadband  signals  were  first 
decomposed  into  33  narrowband  components.  100 
snapshots  were  taken  for  each  of  the  50  runs. 
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An  initial  estimate  o£  18°  was  used.  Two  cases 
were  considered  for  ESPRIT.  In  the  first  case, 
overlapping  subarrays  were  used.  Tims,  the  two 
subarrays  were  composed  of  the  first  and  last 
seven  sensors,  respectively.  In  the  second  case, 
non-overlapping  submrrays  were  used.  Here  4 
pairs  of  sensors  were  generated  by  using  ad¬ 
jacent  sensors.  The  saeple  aean  of  the  estimates 
are  shown  in  figures  1  and  3  for  the  sources  at 
16*  and  22*.  The  corresponding  sample  variances 
are  shown  in  figures  2  and  4  along  with  the 
Craaer-Rao  lower  bound  (CRLB).  In  the  figures, 
the  aoving  window  is  denoted  by  (1),  ESPRIT 
(case  1)  by  2,  and  -ESPRIT  (case  2)  by  (3).  (4) 
represents  the  CRLB.  For  signal  to  noise  ratios 
greater  than  10  dB,  the  three  methods  produce 
comparable  results  for  the  sample  mean.  For  SNR 
lower  than  10  dB,  moving  vindow  appears  to  have 
the  largest  bias.  Hovever,  for  all  value  of  SNR 
between  5  and  30  dB,  the  moving  window  has  the 
smallest  variances.  It  is  concluded  that  the 
moving  window  vorks  well  for  broad  band  sources 
when  applied  in  conjunction  vith  CSS. 


Il, 


Angle= 16° 
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Anlge=20° 


Figure  4 
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